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Abstract. The diamond and completely bounded norms for linear 
maps play an increasingly important role in quantum information sci- 
ence, providing fundamental stabilized distance measures for differences 
of quantum operations. Based on the theory of completely bounded 
maps, we formulate an algorithm to compute the norm of an arbitrary 
linear map. We present an implementation of the algorithm via Maple, 
discuss its efficiency, and consider the case of differences of unitary maps. 



1. Introduction 

The need for physically significant and computable distance measures for 
quantum operations and channels is of fundamental importance in quantum 
information science [20]. Most importantly, it is often necessary to determine 
how far apart two quantum operations, represented by completely positive 
maps, are from each other in some meaningful sense. The diamond norm 
was introduced in [16] for this purpose. It arises from physical considera- 
tions and satisfies the important stability property desired for such measures 
[H [11] . Interestingly, the diamond norm is intimately related to the norm 
of complete boundedness, a notion that has been studied in operator theory 
for different reasons over the past four decades [21]. On finite dimensional 
Hilbert space, every linear map has a finite completely bounded (CB) norm. 
Thus, CB maps are precisely the linear maps in the finite dimensional case. 
In operator theory, CB maps are the natural maps between certain objects 
called operator spaces. Computing the norms of CB maps between cer- 
tain operator spaces introduced in [26] has provided the impetus for recent 
progress on multiplicativity conjectures for quantum channels [9l ll0l[12| l31]. 
CB maps and norms have also arisen in a wide variety of other recent investi- 
gations in quantum information science, including pT4l fT5l [T8l [T9l [25l [271 [30] . 
though the CB terminology has not always been used. 

In this paper, based on the classical and contemporary theory of com- 
pletely bounded maps, we formulate an algorithm to compute completely 
bounded and diamond norms for arbitrary linear maps on finite dimen- 
sional Hilbert space. Along the way, we also provide a brief introduction 
to completely bounded map theory, including the generalized Stinespring 
theorem and Choi-Kraus representation for such maps. We then present an 
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implementation of the algorithm via Maple. We also discuss the algorithm's 
efficiency, and note how it is potentially optimal. Our approach to comput- 
ing these norms is distinct from other known approaches, such as the use of 
semidefinite programming |29] . 

In the next section we recall basic properties of the diamond and com- 
pletely bounded norm, showing how to interpolate between the two. This is 
followed by the introduction to completely bounded maps; our presentation 
here is motivated by that of [21] . In the penultimate section we describe the 
theoretical formulation of the algorithm, and apply it to derive a geomet- 
ric formula for the case of differences of unitary maps. In the final section 
we exhibit code for the Maple implementation of the algorithm, giving an 
explanation of each subroutine. 

2. Linear Maps and Stabilized Norms 

We shall write M n for the set of n x n complex matrices, and regard it as 
the set of operators acting on an n-dimensional Hilbert space represented as 
matrices in a given orthonormal basis. Quantum operations or channels are 
represented by linear maps (j) : M n — » that are completely positive and 
trace preserving (in the Schrodinger picture) or unital (in the Heisenberg 
picture). The dual map <jy : — > M n is defined via the Hilbert-Schmidt 
inner product Tr{cf>(A)B) = Tr(A<ft(B)). 

In quantum information one is often interested in properties of differences 
(j) — if) between pairs of quantum operations. Such a difference is still a 
linear map (a "superoperator"), though not necessarily completely positive. 
In fact, every linear map can be decomposed as a linear combination of 
at most four completely positive maps. This parallels the corresponding 
statement about general operators and positive operators on Hilbert space, 
though the proof is more delicate [321 [21] . 

The 1-norm of a linear map <p is given by \\4>\\\ = supi|_ x -|| 1 <i ||<^P0lli> 
where \\X\\\ = Tr|X|. The operator norm of cj) is \\4>\\ = sup||_x-||<i ||</>(X)||, 
where \ \X\\ = supimi.^ ||-X"£||, and £ ranges over the unit ball of the domain 
Hilbert space for X. Every norm \ \\<p\\\ defines a distance measure d(4>,ip) = 
1 1 \<p — ip\ 1 1. Neither of the distance measures defined by these norms satisfies 
the stabilization property for distance measures of superoperators [16^ [TT, fTTj : 

d(id m (g) tf>, id m g) ip) = d((f), tp) Vm > 1. 

This property implies that the distance between quantum operations is un- 
affected by any ancillary quantum system that is independent of the original 
system. 

The diamond norm is defined in |16[ [T] through partial traces, but it is 
shown to be equivalent to the following quantity. 

Definition 1. For a linear map (f) : M n — ► M^, define 

\\(t>\\o = \\id n ®4>\\i. 
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Though not obvious, the stabilization property is satisfied by || • || . One 
way to see this is through a connection with the completely bounded norm, 
to which we now turn. 

Definition 2. For a linear map 4> ■ M n — > M k , define 

||(^|| cfe = sup ||id m <g> <j)\\, 

m>l 

which we call the completely bounded norm or the cb-norm. 

As a convenience, we adopt the notation <f> m = id m ®4>- It is easily checked 
that ||0 m || < ||0 m +i|| and \\4> m \\ < [21 Chapter 1]. Note also that 

||W0|| C 6 = \\4>\\cb = ||0W|| C 6 for every unitarily implemented map U{-) = 
U(-)U^. Furthermore, we have ||0i <8> 4>2\\cb = H^illcfell^Hcfe) and ||id|| c h = 1. 
The identification in Theorem [3] shows the corresponding properties hold for 
the diamond norm. 

It is possible to relate the completely bounded norm and the diamond 
norm as follows. We first note that by a theorem of Smith [281 [21] the 
cb-norm stabilizes in the sense that for a map (f> : M n — » M k we have that 
H^llcb = \\idk (8> <p\\- Then we make use of the duality relationship [4j given 
by ||0|| = supuxn^i H^WHi, to obtain 

||0|U = \\id k ®<j)\\= sup \\{id k (8) 4t)(X)\h = W^Wo, 

\\X\\l<l 

since $ : M k -> M n . 

Thus, we have the following, which also includes an upper bound [28J. 

Theorem 3. Let eft : M n — > M k , be a linear map, then 

U\\c b = U 1 \\» = U k \\<kUl 

In summary, using the fact that these maps appear as dual pairs, we see 
that for tp : M m — > Mj, ||^||o = ||"0 ||c6 — ^IIV'II an( i th & t the stability of the 
diamond norm |16l Q] is the dual version of Smith's stability for the cb-norm 
|28j . A more refined general upper bound is discussed in the next section. 

3. Completely Bounded Map Primer 

We next give a compressed introduction to completely bounded maps on 
arbitrary operator spaces. As noted previously, in the finite dimensional 
case CB maps on M n are precisely the linear maps (f> : M n — > M k . However, 
the important structural results reviewed in this section are best viewed in 
the more general setting of CB maps as in [21]. In any event, if the reader 
wishes to move directly to the algorithm, this section can be skipped save for 
the structural result Theorem [19] and the norm description of Corollary 1201 

Given a (separable) Hilbert space Ti., we denote the set of (bounded) 
linear operators on H by B(7i). Given operators, ly G B(H),1 < i < 
m, 1 < j < n, we identify the m x n matrix of operators, (Tij), with an 
operator from H (n) = H © . . . H (n copies) to W (m) = H® ...H (m 



-1 



N. JOHNSTON, D. W. KRIBS, V. I. PAULSEN 



copies) by regarding vectors in these spaces as columns and performing 
matrix multiplication. That is, we identify M m>n (B(H)) = B{H {n \n^ m) ). 
This endows M mn (B(7i)) with a norm and this collection of norms on B(7i) 
is often referred to as the set of matrix norms on B(7i). 

Definition 4. Let H be a Hilbert space and let M C B(H) be a subspace. 
Then the inclusion, M m ^ n {M) C M„ hn (B(Tl)) endows this vector space with 
a collection of matrix norms and we call, A4, together with this collection of 
matrix norms on M mn (A4) a (concrete) operator space. When m = n, 
we set M n (M) = M n>n (M). 

Thus, an operator space carries not just an inherited norm structure, but 
these additional matrix norms. 

For basic properties of C*-algebras we point the reader to [HJ. C*-algebras 
are defined abstractly, but every abstract algebra is isomorphic to a concrete 
C*-algebra given by a subalgebra of some B(TC) that is closed under both the 
operator norm (|| • ||) and adjoint (f) operation. If A is any C*-algebra and 
7r : A — > B(TC) is a one-to-one f-homomorphism (and hence an isometry), 
then the collection of norms on M mtn (-K(A)) is independent of the particular 
representation n, and hence, the operator space structure of a C*-algebra 
is independent of the particular (faithful) representation. Hence, each sub- 
space A4 C A is also endowed with a particular collection of matrix norms 
and so we also refer to a subspace of a C*-algebra as an operator space, when 
we wish to emphasize its matrix norm structure. We now give the definition 
of a completely bounded map in the general operator space setting. 

Definition 5. Given a C*-algebra A, an operator space A4 C A, and a 
linear map, <j> : M —* B(H), we define (ft n : M n (M) — >■ M n (B(H)) by 
( t ) n(( a i,j)) = ( ( f ) ( a i,j))- We call (ft completely bounded, if 

\\4>\\cb = SUp 110^11, 
n 

is finite. Here \\4> n \\ = sup {\\(p n (A)\\ : A e M n (M), \\A\\ < l}. 

More generally, any time that A4 and J\f are two spaces, both endowed 
with a family of matrix norms, then one can define the completely bounded 
norm of a map (ft '■ -M. — > N in analogy with the above definition. 

Recalling the upper bound of Theorem [3l for maps whose domain is M n 
and range an arbitrary operator space, a result of Haagerup shows that, in 
general, \\ift\\ c b / HVVnll) no matter how large one takes m [HI p. 114], but we 
do have an upper bound. This result is not explicitly in the literature so we 
provide a proof below, that uses some concepts we will introduce in Section 
4 and, perhaps, illustrates their utility. For now, it is enough to know that 
given a finite dimensional normed space X, there exists a constant a(X) 
called the alpha constant of the space, with the property that 



U\\ cb <a(X)U\\ 
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for any map with domain an operator space that is isometrically isomorphic 
to X as normed spaces. Given two finite dimensional normed spaces X, Y 
of the same dimension one has 



where d(X, Y) denotes the Banach-Mazur distance between the spaces. These 
concepts and results can be found in [22J. 

Theorem 6. Let Ai be an operator space and let <p : M n — > Ai be a linear 
map, then \\<p\\ c b < ra \/™||</>||- 

Proof. Let ||-X"||2 denote the Hilbert-Schmidt norm of a matrix X. Since 
< ll-^lk < V^ll^ll' the Banach-Mazur distance satisfies d(M n ,C n ) < 

\fn. 



Hence, a(M n ) < d(M n ,C n2 )a{C n2 ) < y/n~a(C n2 ). Finally, it is shown 
in [22], that for Euclidean space, a(C m ) < y/rn, from which the result 



We next recall the abstract definition of completely positive maps. 

Definition 8. If A is a unital C*-algebra, then a \-closed subspace S C A 
such that 1 € S, is called an operator system. 

Thus, operator systems are operator spaces and have matrix norms. But 
the additional hypotheses guarantee that if we let A + denote the positive 
elements of the C*-algebra, then S is the span of S + = S n A + , which is 
a cone in S. We also have that M n (S) is the span of M n (S) + = M n (S) n 
M n (A) + . The vector spaces, M n (S) together with the cones M n (S) + is often 
referred to as the matrix ordering on S. 

Definition 9. Given a unital C*-algebra A,S <^ A and a map (ft : S — > 
B(7i), we call (j> completely positive, provided that <p n is positive for all 
n, that is provided that (ojj) G M n (S) + implies that (4>(aij)) € M n {B(7i)) + . 

Definition 10. Given a C*-algebra A and M C A an operator space, we 
set Ai* = {g^ : a G M}, which is another operator space. If <ft '■ M B{TL), 
is a linear map, then we define <fi* : Ai* — > BiTL) by 4>*(b) = cj){W)\ which 
is another linear map. 

The following objects allow one to relate much of the theory of completely 
bounded maps to the more familiar theory of completely positive maps. 

Definition 11. Let A be a unital C*-algebra and let Ai C A be an operator 
space, then we define an operator system Sm C M2(A), by 



a{X) < d{X,Y)a{Y) 



□ 



Combining this result with Theorem [3] we have: 
Corollary 7. Let <fi : M n — > be a linear map, then 



|H| c6 <min{fc,\/^}|H|. 
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Details for the next six results can be found in [21]. We briefly sketch some 
of the ideas of the proofs to help indicate the interplay between completely 
bounded maps and completely positive maps. 

Theorem 12. Let A be a unital C*-algebra, M C A be an operator space 
and let <ft : A4 — > B(TC) be linear. Then \\(f>\\ c b < 1 i/ and only if : Sm :— * 
M2(B(7i)) is completely positive, where 

a\,_f\I n 0(a)\ 
^{tf ^(6)t fxlnj- 

In particular, using this identification of completely contractive (||</>|| c b < 
1) maps with "corners" of unital completely positive maps, one extension 
theorem: 

Theorem 13 (Arveson's Extension Theorem). Let S C A be an operator 
system and let 4> : S —* B(7i) be completely positive. Then there exists a 
completely positive map ip : A — * B(7i) that extends 4>, that is, such that 
ip(a) = 4>(a) for every a £ S. 

quickly yields another: 

Theorem 14 (Wittstock's Extension Theorem). Let Ai C A be an operator 
space and let (f> : M. — * B(7i) be completely bounded. Then there exists a 
completely bounded map ij) : A — > B(TL) that extends (ft and satisfies \\ip\\ c b = 

To obtain the second from the first, one first scales <j) so that \\4>\\ c b = 1> 
then applies Arveson's Theorem to extend $ : Sm — > M2(B(H.)), to ^ : 
M2(A) — > M2(B(H)), and then lets if) be the corresponding (1,2) matrix 
corner of *S>. 

A fundamental result for quantum information is Stinespring's classical 
representation theorem for completely positive maps. 

Theorem 15 (Stinespring's Representation Theorem). Let A be a unital 
C*-algebra and let <j> : A — > B(7i) be a completely positive map, then there 
exists a Hilbert space IC, a bounded operator V : Ti — > KL and a unital f- 
homomorphism, tt : A ^ B(K) such that 4>(a) = V^ir(a)V, for every a G A. 

Note that in Stinespring's theorem, we also have that \\4>\\ c b = 110(1) II = 

\\V^V\\ = \\v\\ 2 . 

This form of the Stinespring Theorem is less common in quantum in- 
formation, but the more standard forms can be readily obtained. Suppose 
4> : M n — * Mk is a completely positive unital map. As every representation 
of M n is unitarily equivalent to a multiple of the identity representation |8j , 
7r can be assumed to be of the form n(a) = Le® a, where Le is the identity 
operator on a suitable dilation Hilbert space E. Further, as is unital we 
have I = 4>{L) = V^V , and hence V is an isometry. Thus we may write 4> as 
4>{a) = V^(Le a)V . The dual of this equation yields the familiar partial 
trace form for a quantum channel, <j> {fl) = T^EiYoV^). 
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In a similar fashion to the extension theorem, Stinespring's Theorem can 
be extended to completely bounded maps. 

Theorem 16 (The Generalized Stinespring Theorem). Let A be a unital 
C*-algebra and let <f> : A — > B(7i) be a completely bounded map, then there 
exists a Hilbert space IC, bounded operators V : 7i — * K, W : TL — * K and 
a unital \-homomorphism, ir : A —* B(JC), such that \\(j)\\cb = ||V||||W|| an d 
4>(a) = V^ir{a)W, for every a G A. 

In the finite dimensional case of a completely bounded map (f> : M n — > 
Mfc, the corresponding canonical forms look like (j){a) = V\Ie <8> a)W and 
0t( a ) = Tv E (WaV^). 

The generalization of Stinespring's theorem to completely bounded maps 
also yields the following "polar form" for completely bounded maps. To 
motivate this result note that for operators, if we set \t\ = VWr, then 

(£1 ^.positive. 

Corollary 17. Let A be a unital C*-algebra and let <f> : A — > B(7i) be 
completely bounded, then there exists completely positive maps, <p\ , <f>2 : A — > 
wiift ||^i(l)|| = ||^(1)|| = II^IU, such that $ : M 2 {A) -» M 2 (B(H)) 
is completely positive, where 

fa b\ (Ma) <t>(b) \ 
*\c d)>-\P{c) ch(d) J ■ 

When \\4>\\cb < 1) then 4>i,4>2 can both be taken to be unital completely 
positive maps. Thus, the above corollary is another example of the meta- 
theorem that completely contractive maps are the "corners" of unital com- 
pletely positive maps. 

The generalized Stinespring theorem, unfortunately, has no good unique- 
ness criteria, unlike the usual Stinespring theorem. The difficulty stems 
from the fact that the two completely positive maps, <f>\, 4>2 are not uniquely 
determined by <f>. Generally, there are many possible extensions of the com- 
pletely positive map $ from the operator system 1S4 to M%(A) and this 
allows for a great deal of non- uniqueness. Ostensibly this follows from the 
fact that the Hahn-Banach Theorem plays a key role in the proof. Thus, 
in particular, the above "polar form" of a completely bounded map is not 
unique. 

Just as one obtains the Choi-Kraus representation of completely positive 
maps from M n to Mk by specializing Stinespring's theorem to these algebras, 
one obtains a similar representation of completely bounded maps. 

Theorem 18 (Choi-Kraus Representation Theorem [5l[T7]). Let cp : M n — ► 

Mfc be completely positive, then there exists matrices, A{ G M n 1 < i < nk, 
such that (j)(X) = J2i A l XA i- 
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Theorem 19 (CB Representation Theorem). Let 4> : M n — ► be a linear 
map. Then there exists matrices, A\ G Mk, n , 1 < i < m, and matrices 
Bi 6 M n< k, 1 < i < m, such that 



(1) 



with 



I = \\4>Hl = II Ei Mil II Ei SjBiH and m < nfc. 



It is important to understand the difference between and the usual 
dual map $ considered in quantum information, so let us dwell on this 
point for a moment with A4 = M n ,B{TL) = Mjt, and (f> : M n — ► 
In terms of Choi-Kraus operation elements, if A G Mk : n,B € M nj fc, and 
: M n -> M fc , is defined by = AXB, then 0* : M„'-> M fc is given by 

0*(X) = B*XA*, while ^ : M fe -► M n is given by $(Y) = BY A, and the 
obvious generalization holds true if <p is given by a sum of such terms. 

Another difference between <fi* and <ft arises when considering the CB 
norm. It is easily checked that \\<f> n \\ = \\(4'*)n\\ an d hence that \\<f>\\ c b = 
||^>*||c6- On the other hand, \\4>\\ c b and ||^|| c 6 can be different. For instance, 
in the case of a completely positive, trace preserving map 4>, the dual $ is 
unital (0t(7) = I), so that ||0t|| c6 = ||^t(/)|| = i, whereas \\<f>\\ cb = \\<f>(I)\\ 
could be larger or smaller. 

We shall call any representation <f)(X) = J2i A; t XBi a generalized Choi- 
Kraus representation. Note that if we have any generalized Choi-Kraus 
representation of (j> then, 



fx ... o\ 



<f>(X) = (A 1 ,...,A r . 



X 



\0 







X J 



fBA 



\B m J 



where the term in the middle represents the m x m block diagonal matrix 
each of whose blocks is X and hence, 



Bi 



. Br, 



EMir /2 iiE^ii 1/2 ' 



which explains the asymmetry in the roles of the A's and B's. 
This also leads to the following very useful result. 

Corollary 20. Let 4> ■ M n — > Mk be a linear map, then 

c 6 = ||0 t ||o=inf{||5:Ml| 1/2 |lE^I| 1/2 }' 



where the infimum is taken over all generalized Choi-Kraus representations 
of*. 
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4. Computation and Estimation of the CB/0 Norm 

In this section we present the theoretical formulation of the algorithm and 
use it to derive a geometric formula for mixtures of pairs of unitary maps. In 
the case of a completely positive map eft, Theorem [T9l coupled with the Choi- 
Kraus representation theorem shows that the CB norm of (ft is exactly the 
operator norm ||</>(/)||. For completeness we provide the direct, elementary 
proof of this fact from [21] with no restriction on the domain of the map. 



Lemma 21. Let P and A be operators on some Hilbert space TL with P 
( P A\ 

positive. Then ( pj > implies that \\A\\ < \\P\\. Furthermore, if P 
is the identity operator then the converse also holds. 

( P A\ 

Proof. To show the forward implication, note that if I D I > then it 



follows that 



P A \ x 
A* Pj \-y 



, A * P s 

* j) > Vx,yG Hs.t. \\x\\ = \\y\\ = 1. 



Thus, (Px\x) + (Py\y) > (Ay\x) + (x\Ay) = 2Re({Ay\x)). Also, the Cauchy- 
Schwarz Inequality tells us that {Px\x) + (Py\y) < \\Px\\ + ||Py|| < 2||P|| 
since ||x|| = ||y|| = 1. Thus, ||P|| > Re ((Ay\x)) Vx, y G TL s.t. ||x|| = ||y|| = 
1. A little thought reveals that this is equivalent to ||P|| > (Ay\x) Vx,y G 
TL s.t. ||x|| = ||y|| = 1, which immediately implies that \\A\\ < \\P\\. 

To show the converse when P = I is the identity operator, we prove by 
contradiction by assuming that \\A\\ < \\I\\ = 1 and that 3x,y G TL such 

that( i j, f\ h 



X ] } < 0. It then follows that ||x|| 2 + ||y|| 2 < 



(Ay\x) + (x\Ay). 

The Cauchy-Schwarz Inequality tells us that (Ay\x)+(x\Ay) < \\Ay\\ \\x\\ + 
\\ x \\\\Ay\\ < 2 ||j4|| ||x|| ||y|| < 2 ||x|| \\y\\. Thus, ||x|| + ||y|| < 2||x||||y||. 
This, however, is impossible since Vx,y G TL it is true that (||x|| + ||y|| ) 2 = 
(||x|| 2 — ||y|| 2 ) 2 + 4||x|| 2 ||?/|| 2 > 4||x|| 2 ||y|| 2 , so ||^|| 2 + ||?/|| 2 > 2 
completing the contradiction. □ 

Theorem 22. Let S C A be an operator system, let B be a C*-algebra, and 
let (ft : S — > B be a completely positive map. Then (ft is completely bounded 
and IMUHMI = ||^(1)||. 

Proof. First note that ||</>(1)|| < \\(ft\\ < ||^|| c b; so we need only show that 

UU< Il0(i)ll- 

Fix n and let A G M n (S) be such that \\A\\ < 1, and let L n be the unit of 
M n (A). Then Lemma [2T] tells us that (^™ j J > 0. Since <ft is completely 

positive, it then follows that fan = (a™/^* f. n ^T\ ) — 0- 

\A l n J \(p n {A) (ft n {l n ) J 

Making use of Lemma I2T1 again shows us that |j</> n (^4)|| < ||(^ n (J n )|| = 
\\4> (1) || . Since this inequality holds for any such A, the proof is complete. □ 
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4.1. The Algorithm. We now turn to the problem of actually computing 
the norm of a linear map 4> : M n — > . By the above results we know that 
to compute the cb-norm we need to do a minimization over all generalized 
Choi-Kraus representations. This turns out to be somewhat more attainable 
than might be imagined and we present an algorithm for computing the 
cb/o-norm of such maps. We first describe the algorithm and then justify 
it later. 

We assume that we are given a map <f> : M n — > some generalized 
Choi-Kraus representation (f>(X) = YliLi MXBi and we wish to compute 

uu = u%. 

Step 1. Find a basis, {Ci, C{\ for the span of {B\, B m } and express 
Bi = Y,di,jCj 

Step 2. Using the expressions for each Bi as a linear combination 
of Cj we may re- write <p(X) = Y?j=iDjXCj. In fact, we have <p(X) = 
>: ; ,l,.Vi>: ;( /,./' i = EjiEidij^XCj. Thus, 



Step 3. Find a basis {E±, E p } for the span of {D±, Di}, express 
each Dj as a linear combination, and repeat Step 2 to obtain 

p 

cj){X) = Y J E i XF i , 
i=l 

where the Fj's are the corresponding linear combinations of the Cj's. 

Remarkably, at this stage it is a theorem that the sets {Ei, ...,E P } and 
{Fi, F p } are linearly independent, and hence this process terminates. 

Step 4. Given an invertible S = (sjj) £ M p with inverse S* _1 = (U j) £ 
M p , let Hi = SijFj, and Gj = ^ UjEi. Then 

||0|| cfe = inf{||^G,GII| 1/2 |lE^H 1/2 }' 

i i 

where the infimum is taken over all invertible matrices S. It is also enough 
to consider positive, invertible matrices for S. 

This algorithm reduces the computation of \\4>\\ c b to a series of matrix 
computations and only the last step might involve a difficult minimization. 



COMPUTING STABILIZED NORMS FOR QUANTUM OPERATIONS 



11 



To begin to justify the algorithm, we begin with the last step. First we 
show that 4>{X) = J^iGiXHi. This can be seen formally, because 



GiXHi — (G\, G p ) 



— (E\, E p )(tijl n ) 



(X ... 
X ... 

\o ... 
fx 

X 





x ) 



(HA 



\0 





(Fi\ 


( s i,jln) 









= Y J E i XF i = 4>{X), 



since the two scalar matrices commute past the direct sum in the middle. 
Note that these scalar matrices behave like "environmental operators" , that 
is, they are operators that act exclusively on the environment of an open 
quantum system. 

Next we need to see that the linear maps from M n to M k , which we denote 
by C(M n , Mk), can be identified with the tensor product, Mk n ®M n k via the 
map that sends an elementary tensor, A®B to the map <f>(X) = AXB. It is 
easily seen that this extends to a linear map, V : M kjn ®M n ^ — > C(M n , M k ), 
that a simple dimension count shows is one-to-one and onto (both spaces 
have dimension n 2 k 2 ). 

We now endow M k ^ n <g> M n ^ with a norm so that T will be an isometry 
when C(M n ,Mj t ) is endowed with the cb-norm. By the CB representation 



M, 



n,ki 



theorem, we see that if we define for U 6 M^, 

: inf{iiEMii 1/2 iiE^n 1/2 }' 



\u\ 



where the infimum is taken over all ways to represent U = Ai <g> Bi as a 
sum of elementary tensors, then we will have that \\U\\h = ||r(i7)|| C 6. 

The above tensor norm is called the Haagerup tensor norm in honor of 
U. Haagerup who was the first to notice the above identification. We write 
Mk, n ®h to denote the tensor product endowed with this norm and 

note that we have just proved that: 

Theorem 23 (Haagerup). The map T : M Kn ® h M n>k -> CB(M n ,M k ) 
defined by F(A (g) B)(X) = AXB is an isometric isomorphism. 

Here we use CB(M n , M&) to denote the space of linear maps from M n to 
Mfc endowed with the completely bounded norm. The above isomorphism 
was greatly extended in work of Haagerup and Effros-Kishimoto to other 
identifications between spaces of completely bounded maps and Haagerup 
tensor products. 
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The above theorem reduces the justification of the above algorithm to 
showing that if <f) = r(?7), then the algorithm correctly computes, ||C/||/i- The 
fact that this algorithm correctly computes \\U\\h for any operator spaces is 
proven in [3]. We outline the key ideas below. For this, we will need a few 
facts about tensor products of vector spaces. 

Recall that if V and W are vector spaces, then every element of V (8> W is 
a finite sum of elementary tensors. The least number of elementary tensors 
that can be used to represent an element u S V (8> W is called the rank of u 
and is denoted by rank (u) . 

Proposition 24. [3] Let u G V (8> W. If u = Y%=i v i ® ^ ^ en P = rank(u) 
if and only if {vi, ...,Vp\ is a linearly independent set and {w\, w p } is a 
linearly independent set. Moreover, if u = X}?=i x « ® 2/» * s another way to 
represent u as a sum of elementary tensors and p = rank(u), then 

span{v\, v p } = span{x%, ...,x p } 

and 

span{wi, ...,w p } = span{y 1 , ...,y p }. 

Proposition 25. [3j Let u <E V ® W. If we apply Step 1 and Step 2 of 
the above algorithm to u = YliLi a i ® h, to obtain u = Y%=i e « ® /*> 
i/ien {ei,...,e p } anrf toZ/ 6e linearly independent sets and hence 

rank(u) = p. 

These facts are easily proved by applying maps of the form / (g) idw and 
idy (8) 5, where / and g are linear functionals to u. 

The remainder of the proof of the justification of the algorithm is to show 
that at each stage, removing the linear dependencies among the elements in 
the sum for u reduces the Haagerup norm. This is best seen at each stage 
of the algorithm. Say at Step 1, when we choose the basis, {C\, C{\ and 
express, Bi = Y^jdijCj, if we first polar decompose the matrix (dij) = 
i w i,j)(Pi,j) where W = (wij) is an m x I partial isometry and P = (pij) 
is an invertible I x I positive matrix, then we have that Bi = ^2iWi,jCj, 
with C{ = Y2j Pi,jC 'j- I n this case, the set {Ci, C{\ is another basis for the 

span of {Ci,...,C,} and EA <5i = EiCjCi. Moreover, using this basis, 
we would obtain another representation for 4>(X) = Y2j=i DjXCj, where 

Dj = WijAi. Again, since P is invertible, the span of {-Di, Df\ is the 
same as the span of {D\, ...D{\. Moreover, since W is a partial isometry, 

one finds that ^2iDiDi < Yli^i-^l- Thus, the infimum of the Haagerup 
norm expression over all linear combinations of the ZVs and Cj,'s which is 
the same as the infimum over all linear combinations of the D^s and Cj's is 
smaller. 

This proves that the quantity defining the Haagerup tensor norm (which is 
the same as the CB norm) must be attained when the coefficients of the gen- 
eralized Choi-Kraus representation are linearly independent, and hence rep- 
resented by some choice of basis for span{E\, ...,E P } and span{F\, ...,F p }. 
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4.2. Example. In quantum information, maps given by the difference of 
two (distinct) unitary maps form the most elementary class of linear, non- 
completely positive maps of interest. We next show how the algorithm can 
be used to derive a simple geometric technique that computes the exact 
stabilized norm for maps in this class. We note this result can be derived 
from a technical result of Herrero ([13], Theorem 3.31), which is proved 
using operator theoretic machinery. Moreover, the result is also stated more 
recently in pQ without proof. Our proof is new and elementary and gives a 
good illustration of the algorithm at work. By the unitary invariance of the 
cb/o norm, observe that we can compute the norm of any map IA — V once 
we know how to compute it for any map of the form IA — id. 

Theorem 26. Let U G M n be a unitary operator on a finite dimensional 
Hilbert space and let <3? : M n — > M n be given by $(X) = UXU^ - X. Then 



cb 



cb 



\&\\ is equal to the diameter of the smallest 



closed disc that contains all of the eigenvalues of U. 

Proof. If U is a scalar multiple of / then $ = so the result immediately 
follows. Thus, we will assume from here on that U is not a scalar multiple 
of /. It then follows from Step 4 of the algorithm that 



1*1 



cb 



inf 



U 



where the infimum is over all invertible 2x2 scalar matrices. 



Now let v = (a c) and w 



U 



2 


faU- 


u\ 




\cU- 


dlj 



(b d y 

2 



so that 



a\ 2 + \b\ 2 + \c\ 2 + \d\ 2 )I - (ab + cd)U ] - (ab + cd)U 

= \\{\\A 2 + ll w H 2 K- 2-Re((v,w)*7)|| 

If we let D = ad — be be the determinant of the matrix, then a similar 
calculation shows that 

2 



D 



-i 



d -b 
—c a 

Thus it follows that 

\m cb =M{\D 



\D 



-2 



|v|| 2 + ||w|| 2 H 



2 J Re((v,w)C7) 



-li 



Ivll 2 + llwll 2 )/ 



2 J Re((v,w)C7)||} 



where the infimum is taken over all 2 x 1 complex vectors v and w. 

Now it is clear that this minimum will be attained when v and w are 
rotated such that minj {i?e((v, w)Aj)} is as large as possible (while keeping 
v and w of fixed length), where Aj ranges over all eigenvalues of U. Thus, 
since multiplying w by e ta will not change \D\, it follows that 

infll^r 1 ! 



cb 



Ivll 2 + ||w|| 2 H- 



2Re({v,w)e ia U)\\} 
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where a is such that the minimum real part of the eigenvalues of e ia U is as 
large as possible. Define r to be this largest minimum real eigenvalue part. 

Now, similar to before, we can multiply w by e l/3 so that | (v, w) | = 
(v, e l ^w), and so it follows that 

||*|U = infdDI^IKIIvll 2 + ||w|| 2 )/ - 2 |(v, w)| rl\\} 

where the infimum is now taken over all 2 x 1 real vectors v and w. It is now 

1 1 2 1 1 1 2 

clear that we can assume without loss of generality that ||v|| + ||w|| =1. 
It also follows from some simple algebra that, given any two vectors v and 
w such that ||v|| / ||w||, the value within this infimum will be made smaller 

1 1 2 1 1 2 1 

by scaling v and w so that ||v|| = ||w|| = ^- 

It then immediately follows from expanding out the terms within the 
infimum that this is equivalent to the following minimization problem 

. . 2\ab + cd\r' 

* u - 



. fl 

= mm < — 



\ad — bc\ 



such that a 2 + c 2 = b 2 + d 2 = ^ 

Now, if r < then it is easy to see that this minimum is equal to 2 by setting 
a = d = and b = c = 0. Thus, it only remains to prove the conjecture in 
the case when r > 0. If r > then it is clear that this minimization problem 
is equivalent to the one we get if we remove the absolute value bars in the 
numerator. 

We now form the Lagrangian of this problem: 

\-2{ab + cd)r , / , 2 A ^ f,2 ,2 1 

If we now set ^ = we arrive at the equation 

(v|w) = ab + cd = 2 (a 2 + c 2 ) (b 2 + d 2 ) r = - 



This, however, implies that 6 = arccos(r), where 6 is the angle between v 
and w. Thus, this problem is minimized by vectors v and w that are each 
of length -4^ and separated by an angle arccos(r). This, however, implies 

that \D\ = ||v|| ||w|| sin# = — r 2 . Plugging this and ab + cd = § into 
the formula to be minimized, we learn that 

\\<S>\\ cb = 2Vl-r 2 



It now is a simple geometric argument that finally shows that this value 
is equal to the diameter of the smallest closed disc enclosing the eigenvalues 
of U, completing the proof. □ 



COMPUTING STABILIZED NORMS FOR QUANTUM OPERATIONS 



15 



5. Implementation of the Algorithm 

The implementation of the algorithm via Maple that we now present is 
split up into several procedures, which will be described as we present their 
code. We have kept variable names within the code as close as possible 
to their counterparts presented in the above theoretical discussion of the 
algorithm. After the code has been presented and briefly explained, we 
discuss its efficiency and provide an example that shows how the code is 
used, comparing the results of the algorithm to known theoretical results. 

The first procedure, RandomPositive, generates a random positive matrix 
with eigenvalues in the interval (evalLower, evalUpper]. This is achieved 
by generating a diagonal matrix with entries contained in that interval and 
conjugating by a random unitary. The random unitary matrix is constructed 
by generating a matrix with random entries from the square with corners at 
and 1 + i and then using the Gram-Schmidt process on its columns. 

> RandomPositive := proc(ndim, evalLower, evalUpper) 

local r,i,RD,RM,V,W,U,P:r:=0: 

RD : =RandomMatrix (ndim , ndim , generator=evalLower+ 

DBL_EPSILDN. . evalUpper , outputoptions= [shape=diagonal] ) : 

while r < ndim do 

RM : =RandomMatr ix (ndim , ndim , generator=0 . . . 1 . 0) + 
I*RandomMatrix (ndim, ndim, generator=0 . . .1.0) : 

r : =round (Rank (RM) ) : 
od: 

for i from 1 to ndim do 

V[i] :=Vector(ndim, (j) -> RM[j,i]): 
od: 

W:=GramSchmidt(convert(V,list) .normalized) : 

U :=Matrix (ndim, ndim, (i,j) -> W[j][i]): 

P[l] :=MatrixMatrixMultiply(MatrixMatrixMultiply(U,RD) , 

Hermit ianTranspose (U) ) : 
P [2] :=MatrixMatrixMultiply(MatrixMatrixMultiply(U, 

Matrixlnverse (RD) ) ,HermitianTranspose (U) ) : 
RETURN (P[l] ,P[2] ) : 
end: 

The procedure IsCPMap determines whether or not the completely bounded 
map that it is given is actually a completely positive map by determining 
whether or not its Choi matrix is positive. This procedure is optional for 
the algorithm; it simply serves to allow the algorithm to compute the com- 
pletely bounded norm of completely positive maps more quickly and more 
accurately than it otherwise could. 

> IsCPMap := proc(CelA,CelB,NumOps,n,k) 

local i , j ,x,Choi ,LgChoi ,MtxUnit : 
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if not (n = k) then 

RETURN (false) : 
else 

for i from 1 to n do 
for j from 1 to n do 

MtxUnit [i] [j] : =OuterProductMatrix(UnitVector (i ,n) , 

UnitVector (j ,n) ) : 
Choi [i] [j] : =add(MatrixMatrixMultiply ( 

MatrixMatrixMultiply(A[x] , MtxUnit [i] [j]),B[x]), 
x=l . .NumOps) : 

od: 
od: 

LgChoi:=Matrix(n~2,n~2, (i, j) -> Choi [floor ( (i-1) /n)+l] 

[floor((j-l)/n)+l] [((i-1) mod n)+l,((j-l) mod n)+l]): 
RETURN (IsDefinite(LgChoi) ) : 
fi: 
end: 

The procedures MakeLinlndep and CellMatricize jointly perform steps 
1 through 3 of the algorithm. They are reasonably straightforward. 

> MakeLinlndep := proc(CelA,CelB, NumOps, n,k) 

local BM,BS,u,x,i, j ,bg,BG,d,C,DOp,CM,CS,v,cg,CG,c,E,F: 
BM:=CellMatricize(CelB, NumOps, n,k) :u:=round(Rank(BM)) : 
BS:=Basis({Row(BM, [1. .NumOps] )}) : 
for x from 1 to NumOps do 
for i from 1 to n*k do 
for j from 1 to u do 

bg[i,j] := BS[j] [i] : 
od: 

bg[i,u+l] := BM[x,i] : 
od: 

BG[x] : =ReducedRowEchelonForm( 

Matrix (n*k,u+l, (i, j) -> bg[i,j])): 

for j from 1 to u do 

d[x,j] :=BG[x] [j,u+l] : 
od: 
od: 

for x from 1 to u do 

C[x] :=Matrix(k,n, (i, j) -> 

add(d[l,x]*CelA[l] [i, j] ,1=1. .NumOps)) : 
D0p[x] :=Matrix(n,k, (i, j) -> BS [x] [j+(i-l)*k] ) : 
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od: 

CM:=CellMatricize(C,u,k,n) : v : =round(Rank(CM) ) : 
CS:=Basis({Row(CM, [1. .u])}) : 
for x from 1 to u do 
for i from 1 to n*k do 
for j from 1 to v do 

cg[i,j] := CS[j] [i] : 
od: 

cg[i,v+l] := CM[x,i] : 
od: 

CG[x] : =ReducedRowEchelonForm( 

Matrix (n*k,v+l, (i, j) -> cg[i,j])): 

for j from 1 to v do 

c[x, j] :=CG[x] [j ,v+l] : 
od: 
od: 

for x from 1 to v do 

E[x] :=Matrix(k,n, (i, j)->CS [x] [j + (i-l)*n] ) : 

F[x] :=Matrix(n,k, (i , j ) ->add(c [1 ,x] *D0p [1] [i,j] ,1=1. .u)) : 

od: 

RETURN (E,F,v) : 
end: 

> CellMatricize := proc(Cel,a,b,c) 
RETURN (Matrix ( a, b* c, (i,j) -> 

Cel[i] [floor((j-l)/c)+l,j-c*floor((j-l)/c)])) : 

end: 



The procedure CBNorm is the main procedure; it calls upon the other 
procedures to compute the CB norm of the given map. The procedure 
begins by determining whether or not the given map is completely positive, 
and if so, returns the map's exact CB norm, using the result of Theorem 1221 
If the map is not completely positive, the algorithm described above begins 
to run. 

Steps 1 through 3 are performed via the MakeLinlndep procedure. The 
minimization in step 4 is approximated by calling upon RandomPositive 
repeatedly to compute random positive matrices with eigenvalues in the 
interval (0, 1] and taking the minimum of the resulting norm estimates. We 
can restrict ourselves to using only positive matrices with eigenvalues in the 
interval (0, 1] rather than all positive matrices because multiplying such a 
matrix by a constant will not change the resulting norm estimate. 
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The inputs to the function are CelA, an array of the map's A opera- 
tors in one of its generalized Choi-Kraus representations; CelB, an array of 
the map's B operators in the same representation; Numlts, the number of 
random matrices to be used to estimate the norm (a higher number will 
produce a more accurate estimate but will take longer to compute); and 
NumOps, the number of A and B operators that there are for the map in its 
given representation. 

> CBNorm := proc (CelA, CelB, Numlts, NumOps) 

local i , x , n , k , CBGuess , NewCBGuess , ST , G , H , GG , HH , GGCB , HHCB , EF : 
n:=ColumnDimension(A[l] ) :k:=RowDimension(A[l] ) : 
CBGuess : =inf inity : 

if IsCPMap (CelA, CelB, NumOps, n,k) then 
CBGuess : =Norm (add (Matr ixMatrixMult iply ( 

MatrixMatrixMultiply (A [x] , IdentityMatrix(k) ) ,B[x]) , 
x=l. .NumOps) ,2) : 

else 

EF : =MakeLinIndep (CelA , CelB , NumOps , n , k) : 

for i from 1 to Numlts do 

ST:=RandomPositive(EF[3] ,0,1) : 

for x from 1 to EF[3] do 
H[x] :=Matrix(n,k, (i, j) -> 

add (ST [1] [x,l]*EF[2] [1] , 1=1 . . EF [3] ) ) : 

G[x] :=Matrix(k,n, (i, j) -> 

add(ST[2] [l,x]*EF[l] [1] , 1=1 . . EF [3] ) ) : 

HH[x] :=MatrixMatrixMultiply( 

Hermit ianTranspose (H [x] ) ,H [x] ) : 
GG[x] :=MatrixMatrixMultiply( 

G[x] , Hermit ianTranspose (G [x] ) ) : 

od: 

HHCB :=simplify (Matrix (k,k, (i, j) -> 

add(HH[l] [i,j] ,1=1. .EF[3]))) : 
GGCB :=simplify (Matrix (k,k, (i, j) -> 

add(GG[l] [i,j] ,1=1. .EF[3]))) : 
NewCBGuess : =Re ( sqrt (evalf (Norm (GGCB , 2) ) * 

evalf (Norm(HHCB , 2) ) ) ) : 
if NewCBGuess < CBGuess then 

CBGuess : =NewCBGuess : 
fi: 
od: 
fi: 
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RETURN (CBGuess) : 
end: 

To use the provided code, run the command with(LinearAlgebra) : in 
Maple and then load in the procedures defined above. As an illustration we 
return to a special case of the class of maps discussed above. 

Example 27. Let U be the 3x3 diagonal unitary matrix with eigenvalues 
e y 4 )^ e y m )^ and e \ 4 ) . Then the following code computes the CB norm of 
the map <f> (X) = WXU - X. 

> Num0ps:=2: 

> Numlts:=100: 

> A:=Array(l . .NumOps) : 

> B : =Array (1 . . NumOps) : 

> A[l] :=DiagonalMatrix( [exp(-5*I*Pi/4) ,exp(-I*Pi) , 
exp(-3*I*Pi/4)]) : 

> A [2] :=IdentityMatrix(3) : 

> B[l] :=DiagonalMatrix( [exp(5*I*Pi/4) ,exp(I*Pi) , 
exp(3*I*Pi/4)]) : 

> B[2] :=-IdentityMatrix(3) : 

> CBNorm (A, B.Numlts, NumOps) ; 

Running this code gives an output of 1.449 in just under 7 seconds. The- 
orem[2E\ tells us however that W^W^ = V% so our algorithm is correct to 
two significant digits. To get a more accurate estimate, we can of course 
increase the number of iterations from 100, and it should be clear how to 
modify this code to find the CB norm of other maps. 

5.1. Efficiency. To look at the efficiency of the algorithm, we assume that 
n = k and consider Steps 1-3 separately from Step 4, as Steps 1-3 need only 
be run once for a given map, while we may wish to run our implementation 
of Step 4 hundreds or thousands of times for the same map. 

First note that the most efficient algorithm that could possibly exist for 
computing the CB norm of a general CB map is 0(n 4 ), which can be seen 
by observing that a general CB map will have about n 2 linearly independent 
generalized Choi-Kraus operators, each of which has n 2 entries that must 
each be read at least once. One can observe, however, that the efficiency 
of Steps 1 - 3 of this algorithm is 0(n 8 ), as we need to apply Gaussian 
Elimination to about n 2 matrices each of dimension about n 2 x n 2 . 

If we know already that the generalized Choi-Kraus operators of our given 
representation are linearly independent, we can simply proceed to Step 4 of 
the algorithm, which has efficiency that can be seen to be 0(n 5 ), as the step 
of computing £V G{G\ and Yli^X^i from the families of matrices {G{\ 
and {Hi} is 0(n 5 ) (if computed using the standard matrix multiplication 
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algorithm) . Since all of the other operations in Step 4 are at least as efficient 
as 0(n 4 ), it follows that if we can find a clever way to compute J2iGiG* 
and ]T\ H*Hi, we can reduce the overall order of this step. 

One of the most obvious ways to improve the order of Step 4 is to use 
an algorithm like the Strassen algorithm or the Coppersmith- Winograd al- 
gorithm [7] to perform our matrix multiplications. Employing these matrix 
multiplication techniques would then reduce the order of Step 4 to about 
0(n 4,807 ) or 0(n 4 ' 376 ), respectively. However, the resolution of either of two 
conjectures in [6] would imply that general matrix multiplication can be car- 
ried out in about 0(n 2 ) time, which would imply that Step 4 can be carried 
out in about 0(n ) time, making our algorithm optimal for maps in which 
we already have a linearly independent representation. 
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